Dynamical Gauge Conditions for the Einstein Evolution Equations 



m 
O 
O 

(N 
Oh! 

< 

00 



> 
O 
(N 



O 
m 
O 

"o ! 

cr 

i ' 



X 



Lee Lindblom and Mark A. Scheel 
Theoretical Astrophysics 130-33, California Institute of Technology, Pasadena, CA 91125 

(Dated: February 7, 2008) 

The Einstein evolution equations have previously been written in a number of symmetric hyper- 
bolic forms when the gauge fields — the densitized lapse and the shift — are taken to be fixed functions 
of the coordinates. Extended systems of evolution equations are constructed here by adding the 
gauge degrees of freedom to the set of dynamical fields, thus forming symmetric hyperbolic systems 
for the combined evolution of the gravitational and the gauge fields. The associated characteris- 
tic speeds can be made causal (i.e. less than or equal to the speed of light) by adjusting 14 free 
parameters in these new systems. And 21 additional free parameters are available, for example to 
optimize the stability of numerical evolutions. The gauge evolution equations in these systems are 
generalizations of the "i^-driver" and 'T-driver" conditions that have been used with some success 
in numerical black hole evolutions. 



I. INTRODUCTION 

The traditional 3+1 approach to the study of the Ein- 
stein evolution equations assumes that spacetime is foli- 
ated by a one-parameter family of spacelike surfaces. The 
spacetime metric is usually decomposed with respect to 
these t = constant surfaces according to 



ds 2 



-N 2 dt 2 + g.jidx 1 + N l dt){dx j + N j dt), (1.1) 



where <?y is the (positive definite) three-metric on the 
t = constant surfaces, and N and N l are called the lapse 
and shift respectively pj. (The x % represent spatial co- 
ordinates on the t — constant surfaces.) If dt is the tan- 
gent vector along which evolutions will be generated, and 
n = d T is the unit normal to the t — constant surfaces, 
then the above definitions imply that d t = Nn + N l di. 
Thus the lapse N measures the rate at which proper time 
r advances (as a function of t) along the unit normals, 
while the shift N % measures the velocity of points with 
fixed spatial coordinates with respect to the unit normals. 

The lapse N and shift N l are therefore descriptions of 
how the coordinates {t, x 1 } are laid out on the spacetime 
manifold, and so in this sense they represent coordinate 
or "gauge" degrees of freedom. The lapse and shift are 
not determined by the Einstein equations, and may be 
chosen quite freely. For example the Einstein evolution 
equations have been written in a variety of symmetric 
hyperbolic forms in which the (densitized) lapse and shift 
can be specified as arbitrary functions of the coordinates 

{*, *<} u a a 1 1 1 1 in m e m 

Since the lapse and shift are not determined by the 
Einstein equations, we have the opportunity and the re- 
sponsibility to specify them in some other manner. We 
may use this freedom in a variety of ways. For exam- 
ple we could use it to simplify the representation of the 
spacetime geometry (as is often done in spacetimes with 
symmetries) [T^, ll3 , to simplify the form of the evo- 
lution equatio ns Il7l to avoid singularities (physical and 
coordinate) p], Il8l Il9| , or to attempt to control the sta- 
bility of numerical evolutions pfl l2ll l23| . 

In this paper we develop systems of evolution equations 
that include the lapse and the shift as dynamical fields. 



These equations together with appropriate versions of 
the Einstein evolution equations form symmetric hyper- 
bolic systems for the combined gravitational and gauge 
fields. Unified hyperbolic systems of equations for the 
evolution of the gravitational and the gauge fields have 
been proposed before. The earliest of these uses harmonic 
gauge conditions that reduce the evolution equations to 
a very simple form 0, 0] > but this gauge has not found 
widespread use in numerical simulations of black hole 
spacetimes. Dynamical gauge conditions have also been 
proposed that convert well-known elliptic gauge condi- 
tions into equations that are hyperbolic when the other 
dynamical fields are considered fixed j^jj, El H^, but 
these equations have never been fully integrated with the 
rest of the Einstein evolution equations to form a unified 
hyperbolic system. Strongly-hyperbolic |26| and more 
recently symmetric-hyperbolic |27| formulations that in- 
clude rather general evolution equations for the lapse 
(but which still keep the shift fixed) have also been pro- 
posed. Here we propose a new symmetric-hyperbolic sys- 
tem that includes dynamical equations for the lapse and 
the shift. Our equations are natural generalizations of 
the "if-driver" and the 'T-driver" equations that have 
been used with some success in evolving black hole space- 
times mmm. 



In Sec.[n]we review the properties of these gauge evo- 
lution equations and in Sec. IIHI combine them with the 
Einstein evolution equations to form a single unified sys- 
tem. In Sec. II VI we show that a 16-parameter family 
of these combined (gauge and Einstein) evolution equa- 
tions is symmetric hyperbolic. In Sec.^we find analyti- 
cal expressions for the characteristic speeds of these new 
systems. These expressions depend on 14 of the 16 free 
parameters. We also demonstrate with specific examples 
that all of these characteristic speeds can be made causal 
(i.e., less than or equal to the speed of light) by making 
suitable choices for the 14 parameters. Finally in Sec. I VII 
we extend the evolution equations by performing a gen- 
eral kinematical transformation on the dynamical fields. 
This transformation depends on 19 additional free pa- 
rameters which leave the characteristic speeds and the 
hyperbolicity conditions unchanged. 
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II. DYNAMICAL GAUGE CONDITIONS 

Our aim is to find equations for the gauge fields that 
allow the spacetime coordinates to adapt dynamically to 
the structure of the evolving spacetime. In particular 
we would like the gauge fields to select coordinates in 
which all the dynamical fields become time independent 
whenever the spacetime itself evolves into an equilibrium 
stationary state. For computational efficiency and ease of 
formulating appropriate boundary conditions, we prefer 
to find hyperbolic rather than elliptic equations for the 
gauge fields. We also prefer hyperbolic equations rather 
than equations of indeterminant type because they have 
a well posed initial value problem. 

The desire to improve the stability and accuracy of nu- 
merical evolutions of Einstein's equations has for many 
years provided the motivation to find intelligent choices 
for the gauge fields P, 0] . Perhaps the most widely stud- 
ied gauge condition of this type is the use of maximal- 
slice foliations for the t = constant surfaces in the 3+1 
decomposition. Maximal slices are defined by the condi- 
tion that the divergence of the normal vector vanishes. 
Maximal slices tend to avoid strong focusing singularities, 
and they allow longer numerical evolutions than do sim- 
pler choices such as N — 1. The mathematical condition 
that a slice be maximal is equivalent to the condition that 
the trace of the extrinsic curvature of the slice vanishes: 
= K = g^Kij. The time evolution of K is determined 
by the standard 3+1 ADM expression 

d t K - N l V. L K = -V^jiV + NK tJ K lj , (2.1) 

where Vi is the covariant derivative compatible with g^. 
Thus the choice of evolving along a foliation of maximal 
slices, each with K = 0, is enforced by imposing an ellip- 
tic equation on the lapse N . This condition for the lapse 
is easily generalized to conditions whose effect is to freeze 
K to its value on an initial surface: = dtK. These a K- 
freezing" conditions also result in elliptic equations for 
the lapse on each time slice: 

= -d t K = V l ViN - NKijK lj ~ N^iK. (2.2) 

The if -freezing conditions have been used numerically 
with some success 28]. One disadvantage is that they 
require the solution of an elliptic equation at each time 
step. This is usually more computationally expensive 
than solving hyperbolic equations, and for the case of ex- 
cised black holes [2(1 123 . 128L l29l I30L l3lL 133 . l33 | it requires 
appropriate boundary conditions |34| to be imposed on 
the excision surfaces. For these reasons, alternatives to 
Eq. I|2.2I) have been studied as well. One possibility is to 
convert the elliptic equation for the lapse into a hyper- 
bolic equation, by adding suitable time derivative terms. 
Thus one might take 

d 2 N + KNd t N = -LiN 2 d t K, (2.3) 

as a gauge condition 01 . The second time derivative 
term df N converts the elliptic equation for N into a hy- 
perbolic equation with characteristic speeds ±^/jEt, while 



the first-order term KNdtN provides dissipation that 
tends to suppress d t N. Gauge conditions of this type 
have been called "AT-driver" conditions [3j| and have 
been used with some success in the numerical evolution of 
black hole spacetimes |2JJ,|2jJ,|23j. A large family of differ- 
ent AT-driver conditions can be constructed from Eq. I|2.3|) 
by adding terms that leave the hyperbolic structure of 
this equation intact. Here we will use as our starting 
point one of these A"-driver equations that admits an 
exact first time integral. Thus we adopt a first-order 
AT-driver condition which can be thought of as the first- 
integral of an equation like Eq. (|2.3II : 

= dtN - N%N + nN 2 + fiN 2 (K - K ), (2.4) 

Here Kq is the arbitrarily prescribed value of K on some 
t = constant surface. Lapse functions that solve this 
equation will also satisfy a damped wave equation that 
is analogous to Eq. 1)2. 3fl . Thus our expectation is that (if 
and) when a spacetime evolves into a time-independent 
state, this choice of lapse will drive the evolution toward 
a slicing in which the trace of the extrinsic curvature K 
takes the time-independent value Kq. 

Next we turn our attention to finding appropriate con- 
ditions for the shift N l . The idea is to use our freedom in 
the shift to select spatial coordinates in which the evolu- 
tion of the spatial metric 9t<7y approaches zero whenever 
the spacetime itself evolves toward a stationary state. 
The time derivative of the spatial metric is given by the 
usual 3+1 ADM expression, 

dtgtj = ViNj + VjiVi - 2NKij = E« . (2.5) 

York pj showed that the integral of the square of Sy + 
Xgijg kl T,ki over a t — constant surface is minimized when- 
ever its divergence vanishes: 

= V, (>:" • A(/-">:) (2.6) 
= V i (V i Ar i + V i A rJ ' + Xg ji V k N k ) 

-2Vj [N(K ji + Xg ji K)] . (2.7) 

This is an elliptic equation for N l whenever A > — 2 [3t| . 
Such a condition selects shift vectors that minimize the 
time derivative of the spatial metric (or more accurately 
the time derivative of the densitized metric g x gij), and 
includes the well-studied minimal distortion shift condi- 
tion (the case A = — §) 0. 

It would be straightforward to convert the shift con- 
ditions of Eq. (|2.7|l to hyperbolic equations by adding 
appropriate time derivative terms, in analogy with the 
derivation of the AT-driver equation for the lapse. How- 
ever, we choose instead to follow a slightly different path. 
Motivated by the work of Alcubierre, et. al 0, |2g we 
consider the quantity 

r=g kl f\ u (2.8) 

where g^ is the conformal metric g^ = g x gij, g = det g^, 
and T l ki is the connection compatable with g^. The 
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quantity T l agrees with the dynamical field used in the 
BSSN 37, 38] formulation of the Einstein equations when 
A = — §. It follows from JUJl that 



■A)/2 ff jij 



(2.9) 



and 



+ 1(1 + A) [E«a,- logg + Sa,^'] . (2.10) 

Thus the T-freezing" condition <9 t (g A :P) = imposes 
an elliptic equation on the shift (for A < 3 in this case). 
This T-freezing differential equation has the same princi- 
pal part as the generalized minimum distortion condition, 
Eq. (|2.6|l . Following Alcubierre, et. al |2lJ we convert 
this elliptic shift condition into a hyperbolic equation by 
adding appropriate time derivative terms, e.g. 



dfN 1 



K Nd t N l = nN 2 d t (g x f l 



(2.11) 



As was the case for the lapse equation, it is possible to 
construct a large family of hyperbolic T-freezing condi- 
tions by adding non-principal terms to Eq. (|2.11(l . By 
adding suitable non-principal terms we can construct 
members of this family that admit exact first integrals. 
So we adopt as our "T-driver" condition one of these ex- 
act first integrals: 







+ KNN l -nN 2 (g x f l 



-tin) 



(2.12) 



Here the time-independent g« Tq is the value of g x T l on 
some particular time slice. Our expectation is that (if 
and) when a spacetime evolves to a stationary state, that 
the T-driver condition will cause the spatial coordinates 
to evolve in a way that tends to minimize the coordinate 
time derivatives of the spatial metric. 

In summary then, we adopt the following K-drrver and 
T-driver conditions for the evolution of the lapse and 
shift: 







= 



d t N - N j d 3 N + fi L N 2 (K - K ) 

+n L N 2 - e L Nd t N\ (2.13) 
d t N* - N>d,jN l - ^ s N 2 {g x r - g x T l ) 

+n s NN l - e s Ng ij djN. (2.14) 



III. UNIFIED EVOLUTION SYSTEM 

The if-driver and T-driver Eqs. I|2.13|l and (|2.14|l were 
each constructed to be first-order hyperbolic equations. 
However, these equations are manifestly hyperbolic only 
when the other dynamical fields {e.g., gij, Kij) are fixed, 
whereas the situation of interest to us is when all fields 
evolve together. So our aim now is to construct a unified 
system of evolution equations for both the gauge and the 
gravitational fields such that the entire system is sym- 
metric hyperbolic. 

The first step is to examine the highest derivative cou- 
pling of the (densitized) lapse and shift to the Einstein 
evolution equations. We use a general form of the equa- 
tions written in the notation of Kidder-Scheel-Teukolsky 
(KST) 12]. These are first-order evolution equations for 
the spatial metric gij, the extrinsic curvature and 
the spatial derivatives of the metric Dkij = ^dkgij- At 
this point we need consider only the highest derivative 
(or principal) parts of the equations: 



d t9ij ~ N n d n9ij + 2g n(i d j) N n , 
d t Kij ~ N n d n Kij - NQtdjQ - N 



(3.1) 

(1 + 2a)g cd 5 n {l 8» j) 



-(i + 09 nd s b (l 6 c j} -(i-09 bc s n it s d 3) 



d n D 



bed i 



S n k S\5 c 



(3.2) 



+g nb 5 c i 6 d j + ■>-</" V'";/,, 
d t Dui - N n d n D kij + g aii d j} d k N a - N 

-hrig nb 9k(iS c j) - hxg nb 9ijS c k + hvg^gmS"^ 
+ixg bc g^s n k ]d n K bc , (3.3) 

where ~ denotes equality of the principal part of the 
equation, and Q = \og(N / 'g a ) is the densitized lapse. 
The parameter a that appears in these equations is part 
of the definition of the densitized lapse Q, while 7, 77, \i 
and £ were introduced by adding multiples of the con- 
straints to the evolution equations (see KST [l^l. 

The Einstein evolution Eqs. I|3.2|l and (|3.3|l couple to 
the second spatial derivatives of the densitized lapse and 
shift. Thus in order to construct a first-order unified 
system, we need to promote the spatial derivatives of the 
gauge fields to the status of independent dynamical fields; 
so let 



(3.4) 
(3.5) 



These conditions are just the X-driver and T-driver con- 
ditions of Eqs. H2.4|) and l|2.12|l except for the addition 
of coupling terms between the equations that are pro- 
portional to ei, and eg. These coupling terms will give 
us more flexibility later in constructing a unified system 
of fully hyperbolic equations for the evolution of all the 
gravitational and gauge fields. For maximum flexibility, 
at this stage we take the 7 parameters A, a*s, kl, k s, 
cl and eg to be completely free and undetermined. 



Using these definitions we express the gauge evolution 
equations Ij2.13|l and (|2.14|l in terms of these new fields. 
Furthermore, we obtain evolution equations for Tj and 
Mk by taking spatial gradients of Eqs. (|2.13|l and l|2.14|l . 
The principal parts of the resulting equations are then 
given by 



dtQ 
dtN' 



(3.6) 
(3.7) 
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d t T t ~ N k d k T i +N{2a- l x L )d i K 

+ N(e L -2a)diM/, (3.8) 
dtM/ ~ .Y A «V.U,' • 2N^ s g im g kl d j D klm 



+e s Ng ik d j T k . 



(3.9) 



In deriving the last two equations we made use of the 
constraints 



■^3 



C\, e 20 [i T jl =0, (3.10) 
2iV- 1 a [ „(7VM fe] l ) =0, (3.11) 



in order to write all the terms involving N k as advection 
terms |45j. 

The system of Eqs. (EHJ-ESJl and (E2>h<ElI> consti- 
tutes a unified system of first-order evolution equations 
for the full set of dynamical fields {gij, Kij, Dkij, Q, N l , 
Ti, M k 1 } as desired. However, this system is not unique. 
We are free to add multiples of the various constraints 
to these equations, thus producing other systems whose 
constraint-satisfying solutions are identical. Motivated 
by the fact that the addition of such constraint terms 
improves the mathematical character of the Einstein evo- 
lution equations 01 > we now additional multiples of 
the constraints to our unified system of equations. In 
particular we modify Eqs. (|3.3|l . (|3.8|l . and (|3.9|) as fol- 
lows: 



d t T = 
dtMS = 



(3.12) 



+ ^ 5 Ng lb g ca C ajbc + ^ 6 Ng ia g oc C abcj 



1 r pCri 
Q ^ajbc 



dtD 



klj 



-^i>wNg k (iCj) a a . 



l jifaNgijCka a 



(3.13) 



(3.14) 



Here the . . . denote the terms in the unmodified equa- 
tions. The new terms (each proportional to a new con- 
stant ipA) include multiples of the new constraints, Cy 
and C nk of Eqs. (|3.10|) and (|3.11|) . as well as multiples 
of the standard Hamiltonian and momentum constraints 
C and Ci, and the constraint C k uj from the fixed-gauge 
Einstein evolution system. These latter constraints are 
defined by Hjj 



C = i [^R — KijK ij + K 2 ] . 



C, = VjK j i- 
Ckiij — 2d[ k Di]i 



(3.15) 
(3.16) 
(3.17) 



Adding constraints in this way is essential for obtaining a 
hyperbolic system of evolution equations. Note that the 
constraints added here are in addition to the constraints 
already included in Eqs. (|3.2I) and (|3.3|) . 

The full unified system of evolution equations, includ- 
ing these new constraint terms, can now be written as 
follows (showing here only the principal parts): 



d t gij ^ N k d k gij, 



(3.18) 



d t Q 
d t N l 
d t Ti 



(3.19) 
(3.20) 



(3.21) 



d t D k 



t-LJkij 



0. 

o, 

N k d k T, + N(2a - n^diK 
+N(e L - 2a)a i M j i 
+^ 1 Nd [k K i] k + ^8^^, 
N k d k M/ + 2N f i s g im g kl d j D klm 
+N[2e s a - fx s (l + \)]g lm g kl d.D^ 
+e s Ng tk d 3 T k + ^Ng tk d {l T k] 

+N[^ n 5fg bc + iP 5 g^ b g^ n 6f 

+Ag t[n g a]{b sf + ^g a[b g n]c S<\d n D abc , (3.22) 
N n d n D kij +N[g b[i S^6i + ^g b(i S§8 a k ] 



d n M n 



-N 



5 n k 5 b i S c j ~ hrig nb g k(l 5 c 3 ) - hxg nh g^ k 



+tV9 bc 9k<tS n i) + hX9 hc 9^5 n k ]d n K bc . (3.23) 



dtKu 



N n d n K l3 - Nd^j) - N 
-(l + C)g nd S b {i S c j) -(l-Og bc S n {i S d j) 
+g nb 5 c iS d j + 2 1 g n ^ c g ij \d n D 



(l + 2a)g cd S n (i 6 b j) 
i 

/,,,/. (3.24) 



These equations constitute a first-order system of evolu- 
tion equations for the dynamical fields {gij, K^, D k ij, Q, 
Ti, N l , M k 1 }. This system depends on 22 freely specifi- 
able parameters: 20 of these parameters affect the prin- 
cipal parts of the equations {a, 7, 77, x, (, tpi, . . . , ipio, 
X, fi jj , fis, ^l, £s}j while 2 additional parameters {kl, 
Kg} are dissipation terms in the gauge equations that do 
not affect the principal parts. We will constrain some of 
these parameters in the following section to ensure that 
the system of equations is symmetric hyperbolic. The 
remaining parameters will be freely specifiable and avail- 
able for other purposes, such as simplifying the resulting 
equations or optimizing the stability of numerical space- 
time evolutions. 

These evolution equations for the dynamical fields {gij, 
, D ki j ,Q,Ti,N l , M k 1 } also imply evolution equations 
for the various constraints of the system. In Appendix 1X1 
we derive these constraint evolution equations, and show 
that if the constraints are exactly satisfied initially then 
they will continue to be satisfied as the system evolves. 



IV. SYMMETRIC HYPERBOLICITY 

The unified system of evolution equations H3.18H3.2l|) 
derived in Sec. IHIIcan be written in the form 



d t u a + A ka pd k u t3 ~0 



(4.1) 



where u a is the collection of dynamical fields: u a = {gij, 
Kij, D k ij, Q, Ti, N l , M k {. A first-order system such 
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as this is called symmetric hyperbolic if there exists a 
symmetric positive-definite "symmetrizcr" S a p on the 
space of dynamical fields such that A h a g = S ai _ l A k ^^ is 
symmetric for all k: A k a ^ = Ap a . Symmetric hyper- 
bolic systems [3!| have well-posed initial value problems, 
real characteristic speeds, complete sets of characteris- 
tic eigenvectors, and other nice mathematical proper- 
ties such as the existence of associated canonical energy 
norms. 

We now explore the conditions under which the uni- 
fied evolution equations of Sec. Illll arc in fact symmetric 
hyperbolic. We assume that the symmetrizer S a p can be 
written as a function of the metric gij and various con- 
stant parameters. In particular we consider the follow- 
ing general symmetrizer which we express as a quadratic 
form |47| . 

dS 2 = S af3 du a du p 

= A 1 dG 2 + A 2 g lk g 3l d~ 9lJ dg kl + A 3 dQ 2 
+A i9lj dN l dN J + B 1 dK 2 + B 2 g ik g jl dKij dK u 
+C ig kl g la g^dD {kl]} db (lab] 
+C 2 g kl g la g^ [dD klJ - dD (klj] ] [dD lab - dD (lab) ] 
+C 3 g ij dD\ dD) + C 4 g i3 dD 2 dD 2 + 2C 5 g i3 dD}dD 2 
+E l g lj dT l dTj + 2D 1 g lJ dT l dD) + 2D 2 g lj dT, dD 2 
+E 2 dM 2 + ±E 3 [g l3 g kl + 8\5 k ] dM k l dW 
+iE 4 [g l0 g kl - 5\8 k ] dM k * dW 
+2D 3 dMdK + 2D i g ik 5 J l dk ij dM k l . (4.2) 

Here dG, dK and dM are the traces of dgij, dKij and 
dM k respectively, and dg^ , dK^ and dM k l are their 
trace- free parts. The two traces of dD k ij are defined by 

dD\ = gi k dD ljk (4.3) 
dD 2 = g jk dD kij , (4.4) 

and its trace-free part, dD ki j, is 

dD kij ee dD klJ + \[dD\ i g j ) k -2dD 1 k g ij 

+dD 2 g l3 -'idD 2 t g 0)k \. (4.5) 

The quadratic form 14. 2|) is positive definite iff the 
parameters corresponding to diagonal symmetrizer ele- 
ments {Ax, A 2 , A 3 , A 4 , Bx, B 2 , C u C 2 , C 3 , C 4 , E u E 2 , 
E 3 , E 4 } are positive, and certain inequalities are satisfied 
by the parameters corresponding to the off-diagonal sym- 
metrizer elements {C5, Dx, D 2 , D 3 , D4}. Some of these 
off-diagonal inequalities are simple, i.e. D 2 < E 2 B\, and 
D 2 < B 2 E 3 . But the inequalities involving the other off- 
diagonal parameters {C5, Dx, D 2 } are less transparent. 
The needed condition is that the 3x3 matrix 

C 3 C 5 Dx\ 

C 5 C 4 D 2 (4.6) 
Dx D 2 Ex J 



is positive definite. The most straightforward way to en- 
force this condition is to use the fact that a matrix is pos- 
itive definite iff it admits a Cholesky decomposition [4(| . 
By writing the Cholesky decomposition of Eq. (|4.6[1 in 
terms of new parameters Fa, we obtain 



c 3 


= F 2 + F 2 +F 2 , 


(4.7) 




= F 2 + F 2 , 


(4.8) 


c 5 


= F 2 F 4 + F 3 F§, 


(4.9) 


Dx 


= F3F6, 


(4.10) 


D 2 


= F 5 F e , 


(4.11) 


Ex 


- Fl 


(4.12) 



Given these expressions for {C 3 , C4, C5, Dx, D 2 , Ex}, 
the matrix in Eq. (|4.6[) is positive definite for arbitrary 
F A , so long as Fx ^ 0, F 4 ^ 0, and F 6 ^ 0. 

It is straightforward (but tedious) now to evaluate the 
conditions on the various parameters needed to guarantee 
that the matrices A k a p = Sa^A^ are symmetric in a 
and [3 for all k. After lengthy algebraic manipulations 
we find that the following conditions are necessary and 
sufficient to guarantee that the A k a fj are symmetric: 

B 2 = (Cx+2C 2 )/3, (4.13) 
C = -Cx(E 3 +D 4 )/[E 3 B 2 -D 2 ], (4.14) 

a = [3(3D3 + 2Ai) [-3(2C*i + 5C 3 + 5C 4 + 10C* 5 )£i 
+5(£>i + D 2 ){3Dx + 3D 2 + 3D 3 + 2D 4 )} 
-(3E 2 + 2E 3 ) [9Ex(2Cx + 5C* 3 + 5C 4 + 10C 5 ) 
+5(Dx + D 2 )(9{Bx -Dx- D 2 ) + 2C X + 4C 2 )] } 

/[-10Ex{3E 2 + 2E 3 ){9Bx + 2C X + 4C 2 ) 
+30Ex{3D 3 + 2D 4 ) 2 ], (4.15) 

/is = (I6C1 +5C*2 + 15C 4 )/30(£;3 + -B4), (4.16) 
^ 8 = 2(Cx-C 2 + 3D 4 )/3C 2 , (4.17) 

e s = {5Ex{3E 2 + 2E 3 )(l- X)ll s -6ExCx 

-15(C 3 + C* 4 + 2C 5 )Ex - lQEx{3D 3 + 2D 4 )a 

+5(Dx + D 2 ) [3(Dx +D 2 + D 3 ) + 2D 4 ] } 

/b{3E 2 + 2E 3 ){Dx+D 2 -2aEx), (4.18) 
li L = 2a + [3Bx + 2B 2 - 3(Dx + D 2 ) 

-(3D 3 + 2Dx)e s ]/3Ex, (4.19) 
e L = 2a - [3{Dx + D 2 + D 3 ) + 2D 4 

-{3E 2 +2E 3 )e s ]/3Ex, (4.20) 

ijj 3 = -(l + E 3 /E 4 )e s + (D 2 + D 4 )/E 4 , (4.21) 
V> 4 = {2Ci-5C 2 -15[C5 + (1 + 2(t) j D 4 ] 

-15(^3 + E 4 ) [(1 + X)fi s - 2ae s ] }/WE 4 , (4.22) 

^5 = (Cx+2C 2 +3D 4 )/3E 4 + (Cx-(D 4 )/E 3 , (4.23) 

fa = (Cx +2C 2 + 3D 4 )/3E 4 - (Cx - (D 4 )/E 3 , (4.24) 

7 = {[3Ci-15(C 3 + C 5 ) + 30aDx](D 3 + E 2 ) 
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+5(65! + B 2 )E 2 - 5(6D 3 + D 4 + 3e L D 1 )D 3 
-15/iiZ>i£ 2 + 10(D 4 E 2 - D 3 E 3 )(l + A)ms 
-10<t(3D 3 + 2D 4 )D 3 + 20ae s (A3£ 3 - D 4 £ 2 ) 
+10(7(3^ + 2B 2 )E 2 \/45(Dl - B X E 2 ), (4.25) 

X = [[2{C 3 + 2C 5 )E 1 -2{D l + 2D 2 )D 1 ] 

x [30//sZ3 4 -E'i + (16Ci + 5C* 2 + lbC i )E l 
+ 15(5 2 -D 2 - e S Di)D 2 ) 
+ [(2C 4 + C 5 )Bi - (Di + 2D 2 )D 2 ] 
x [30(£> 2 + e S D i )D 1 - 2(8d - 5C 2 + 15C 6 )J5i 
+ 15Z? 4 Si(2^s + 2A/i5 - ^5 + V»6 - 4ae s ) 
-30S 2 (£>i - 2(7^ + CEi)] }/75^i [CsD 2 
+C 2 £i - 2C 5 D 1 D 2 + (5? - C 3j Bi)C 4 ] , (4.26) 
T] = [15(2 - X )D 2 2 + 15(265^4-3x^1)^2 

+ 10C 2 ^i + 15(V> 5 - V* - ^s)D i E 1 
+ 15(x - 2)C 4 Sx - 30(D 2 + Ei - CEi)B 2 
+45 X C 5 £i]/15[(D 1 + 2D 2 )D 2 - (2C 4 + C 5 )£i] , 

(4.27) 

V^2 = {-16<7i(C 4 Di + 2^2-2^1-^2) 

-5C 2 [C 4 5! + C 5 (Dj - D 2 ) - C 3 D 2 ] (4 + 3^ 8 ) 
+30(C 3 C 4 - Cl) [Di + (E 4 - E 3 )e s + E^ 3 ] 
+60(C 3 D 2 - C 5 Di)(E 3 - Ei)fi s 

-(C 4 D! - C 5 D 2 )t}/30[£i(C 5 2 - C 3 C 4 ) 

+C 3 D2 + C 4 D 2 -2C 5 D 1 D 2 ], (4.28) 
^ 7 = [2Ci - 5C 2 + 15C 4 + 30C 5 - 60D 3 - 9O7D3 
-10D 4 - 10(E 3 + 6£ 2 - 3Ei)[i s 
+30(e L -2<j)D 2 ]/45E 2 , (4.29) 

ih = [-2B 2 -(S X + v)Di 

+ (2- X -2 V )D 2 + 2Die s ]/E 1 , (4.30) 

V'io = {(3C 3 + C 5 )[60(£ 4 -£ 3 )// s -5C 2 (4 + 3V>8) 
+32Ci + 30D 2 ^ 2 - 15C 4 V 8 + 30C 5 (2 + 
-(C 4 + 3C 5 )[l6d + 5C* 2 (4 + 3^ 8 ) + 30D X V>2 
-15C 5 i> 8 + 30C 3 (2 + V 8 )+T]} 

/150(C 3 C 4 - Cf ), (4.31) 
^ 9 = [32Ci- 5(4 + 3^ 8 )C 2 - 15(^8 + 4^10)^4 
+30(2 + Vs - V>io)C s + 60(£ 4 - E 3 )n s 

+30^2D 2 ]/30(C 4 + 3C 5 ), (4.32) 

where T is given by 

T = 30(1 + A)(£ 3 - £4)^3 + 30D 4 (2cr-C) 
-15E 4 (2i/j 4 + + i>6 ~ 4aes) 
+15£ 3 (V>6 - ^5 - 4ae s ). (4.33) 
These conditions determine the 20 parameters {B 2 , ipi, 



• ■ ■ , "010, cr, 7, 7?, x, C, Ml, MS, e L , e s } in terms of the 15 
parameters {A, By, C\ , C 2 , C 3 , C 4 , C 5 , D 4 , D 2 , D 3 , D 4 , 
£2, £3, S 4 } jlq- Writing the conditions for sym- 
metric hyperbolicity like this is a particularly convenient 
way to parameterize these evolution systems. The pa- 
rameters {A, Bi, Ci, C 2 , C 3 , C 4 , C 5 , Di, D 2 , D 3 , D 4) 
i?i, i?2, E 3 , E 4 } can be chosen freely except for the sim- 
ple inequalities needed to guarantee the positivity of S a p- 
We note that the evolution system is invariant under an 
overall scaling of the symmetrizer. Thus without loss of 
generality we will set C% = 1, so there are really only 
14 freely specifiable parameters that affect the evolution 
equations. We also point out the following nice feature 
of this way of parameterizing these equations: By us- 
ing the symmetrization conditions in Eqs. (|4.13|l " (|4.33|) 
to determine the parameters that actually appear in the 
evolution equations {ipi, . . . , ip 10 , cr, 7, r), x, Ml, Ms, 
£l, £5} we are guaranteed to have a system that has only 
real characteristic speeds, a complete set of eigenvectors, 
etc. This same parameterization technique has been used 
by Frittelli and Reula 0, and can also be used to pro- 
vide a more convenient and complete characterization of 
the symmetric hyperbolic subset of the original fixed- 
gauge KST equations. We summarize this approach to 
the KST equations in the Appendix. Finally we note that 
the four symmetrizer parameters {A\, A 2 , A 3 , Ai} do not 
enter any of the symmetry conditions. So while these pa- 
rameters can be chosen quite freely (ensuring only that 
they are positive) they do not seem to play any important 
role in determining the dynamics of the system. 



V. CHARACTERISTIC SPEEDS 

The evolution equations for the full system of fields — 
including the gauge fields — have been put in a first-order 
form in Sec. 11111 The characteristic speeds in the di- 
rection are defined as the eigenvalues of the matrix 
£, k A ka p that appears in Eq. (|4.1|) . The unit one-form 
specifies the direction of propagation. The charac- 
teristic speeds associated with the fields {gij,Q, N 1 } are 
very simple. In the frame of the hypersurface-normal 
observers, the characteristic speed associated with the 
propagation of g^ is v — 0, while the speed associ- 
ated with the propagation of the gauge fields {Q, N 1 } 
i sv = -£, k N k /N gg|. 

In order to evaluate the characteristic speeds as- 
sociated with the other dynamical fields of this sys- 
tem it is convenient to transform to an irreducable 
representation of the space of fields. In this ba- 
sis the matrix £, k A ka p becomes block diagonal and 
hence its eigenvalues become much easier to evaluate. 
The irreducable representation of the remaining dy- 
namical fields {Tj, M k l , Kij, Dkij} consists of pro- 
jecting them onto the scalars {T£\ M fc l £ fc £;, M k i P k i , 
Kigtf, A',, /"*'. D kij e?t. j , D kij t*P ij , D Mj P k *?} 
(where Pij — gij — is the projection ten- 

sor onto the two-space orthogonal to the trans- 
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verse vectors {T.PV,, M k ^ k P^, M k l ^P k ° , l\,iX'l' kj - 
Dku^eP 13 , DmP kj ^\ D kll P k( - l P l », D kll P k iP a }, 
the symmetric transverse traceless tensors {M k iP kl ij, 
K kl P k \ 3 , D mk ^ m P a l3 , D mk iti l P mk l3 } (where P kl t] = 
P k iP l j — ±P kl Pij), the antisymmetric transverse tensors 
{AI k l P k [ i P^i, K k iP k [iP l jA, and finally the transverse 
traceless part of D k ij ■ 

The scalar parts of the dynamical fields form an eight 
dimensional subspace, and this 8x8 block of £, k A ka p 
decouples from the others. This block depends on the 
dynamical KST parameters and the other parameters in- 
troduced in Sec. IIIII that describe the dynamics of the 
gauge fields. We find that the eight characteristic speeds 
(relative to the normals of the hypersurface) can be rep- 
resented as, 

uf 1± = A S1 ±B S1 , (5.1) 
v% 2± = A S2 ±B S2 , (5.2) 

where 

An = i [ml + (1 - + e s e L ] , (5.3) 
B%i = A 2 sl + (l-X)(e L -fi L )fis, (5.4) 
A S2 = i[{l + 2rf)(2 + 2 X -v)-vC] 

+ ^(2 + ^ + 2^10)^5 -V'e) 

+ i^ 7 (2 + 3^ 8 -4^ 9 + 2^io), (5.5) 
B% 2 = A 2 S2 -i[(l + x )(2 + i> s + 2ip w )+r)(i> s -2i> 9 )] 

x[(l + 2 7 )(^ 5 -V 6 )-2C^ 7 ]- (5.6) 

The transverse vector parts of the dynamical fields con- 
stitute two identical eight dimensional subspaccs. The 
eight characteristic speeds for each of these blocks are, 

v 2 V0 = 0, (5.7) 
Vvi = Ms, (5-8) 
v 2 V2 ± = Av2±B V 2, (5.9) 

where 

A V 2 = +^3) + #4(2 + ^8-3^9-^10) 

+|^5(2 + to - ^9) - ^6(2V>9 + 3^io) 
+ ^r;(l-3C-4a)-ix(l + 6a), (5.10) 

By 2 = Ay 2 + xA x + t]A n + A, (5.11) 
A x = ^ (2 + V 8 )[(1 + 3C)V4 + (l + 4£T + CMs 
-(2<t-C)^ 6 ]-^io[(5-9C)V'4 
+ (1 - 4a - 3C)Vb + (4 + 14ct - 3C)to] 

-^[(1 + 6(7)^+3^4 + ^5+^8], (5.12) 

A„ = [(l-4<7 -30^3 -2^4-3^ 6 ] 

+£(2 + to) [3(C - 2a)to - (1 - 4a - 3C)to 
+ (5C - l)to] + 3^9 [(1 - 4(7 - 3C)Vs 
+ (5-9C)^ 4 + (4 + 14a-3C)to], (5.13) 
A = 3 L^iV'io[(l-4a-3C)to-2^4-3to] 
+ ^Vi(2 + to) [^5 + (2d - c)to + to] 
-^to^9 [^5 + (1 + 6(7)to + 3^4 + to] • (5.14) 



We note that the parameter /is introduced in Eq. (|2.14() 
represents one of the characteristic speeds of this system, 
as expected. 

The characteristic speeds of the two identical 4- 
dimensional spaces of symmetric transverse traceless sec- 
ond rank tensors are: 

v \tt)\ = 1; (5.15) 
«(TT)2 = |(^5-^ 6 )(2 + to)- (5-16) 

The characteristic speeds of the 2-dimensional space of 
antisymmetric transverse-traceless second-rank tensors 
are: 

V[TT] = |(^5+^ 6 )(2 + to)- (5.17) 

And finally the subspace consisting of the transverse 
traceless part of D k ij has only one characteristic speed, 
and this vanishes. 

These expressions determine the characteristic speeds 
in terms of the parameters {to? "02, ■■•} that define the 
form of the evolution equations. The speeds can also 
be re-expressed in terms of the symmetrizer parameters 
through Eqs. (|4.13l) - (|4.33|) . The characteristic speeds are 
therefore functions of the 14 parameters {Bi, C 2 , C3, C4, 
C 5 , Dt, D 2 , D 3 , D 4 , Ei, E 2 . E3, Ei, A} that can be 
specified (almost) freely as discussed in Sec. IIVI 

Although most of the characteristic speeds depend on 
the parameters that define the system of evolution equa- 
tions, several of the speeds are independent of them. For 
instance, twelve eigenvectors have characteristic speed 
zero, and four eigenvectors have characteristic speed ±1 
(the speed of light in our units); the former correspond 
to gauge-dependent fields, while the latter must be the 
incoming and outgoing fields corresponding to the two 
physical gravitational degrees of freedom. The remain- 
ing speeds, the adjustable ones, must correspond to var- 
ious gauge-dependent and therefore basically unphysical 
characteristic fields. 

In the past it has been considered most natural 0, 0, 
l+H to set any adjustable speeds in the Einstein evolution 
equations to one (the speed of light) or zero with respect 
to the t = constant surface normals. Our experience, 
however, is that the instabilities limiting evolutions of 
black hole spacetimes (with excision) often occur in out- 
going characteristic fields that propagate at the speed of 
light just outside the event horizon [13|. Excitations in 
such fields remain in the computational domain for long 
periods of time and therefore have the opportunity to 
grow large. It therefore might be better to set the ad- 
justable characteristic speeds to values significantly less 
than the speed of light for evolutions of black hole space- 
times. 

We have not been able to show that the adjustable 
characteristic speeds can be set to arbitrary values by 
adjusting the available parameters, and in fact it appears 
likely that this is not possible. In particular we have 
not been able to find parameter values that make all of 
these speeds equal to unity or zero. However, we have 
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shown that parameter values can be chosen to make all 
of the adjustable characteristic speeds causal (i.e. less 
than or equal to the speed of light) |5(ij . To provide a 
specific example, we have have found parameter values 
that make the characteristic speeds take the following 
"simple" values, 



= 



J (TT)2 



J S1- 



J S1- 



J [TT] > 



,J S2- 



J S2+ 



V V2-1 



V 



V2+i 



(5.18) 
(5.19) 
(5.20) 



to any desired accuracy. The approximate values of 
the symmetrizer parameters needed to achieve these 
characteristic speeds are: B\ = 7.17, C\ = 1.00, C 2 = 
2.57, C 3 = 8.68, C 4 = 3.95, C 5 = -3.81, D 1 = 5.36, D 2 = 
4.86, D 3 = -10.78, D 4 = -2.04, E x = 44.64, E 2 = 
19.39, E 3 = 3.65, E± = 2.22, A = -0.33. These sym- 
metrizer parameters also determine the parameters that 
define the explicit form of the evolution equations; for 
this example the latter parameters have the follow- 
ing approximate values: Vi = 0.13, = —0.34, tfj 3 = 
4.54, i/j 4 = -0.92, V>5 = 0.00, = 0.00, ip 7 = -0.90, Vs = 
-2.00, ^9 = -0.23, ip w = 0.27, 7 = -0.76, a = 0.50, C = 
-0.49, rj = 0.93, x = -0.43, fi s = 0.50, fi L = 0.63, e s = 
— 1.24, 6l = 0.44. We see that all of the characteristic 
speeds in this example are causal, and the various pa- 
rameters that determine the evolution equations are all 
of order unity. In another example, we explored the pos- 
sibility of making all of the characteristic speeds which 
appear in Eq. (|5.18|) (|5.20l) as small as possible. We 
found that it was only possible to make the squares of 
all these characteristic speeds smaller than about 0.29. 
Thus it is relatively easy to find examples of these evolu- 
tion equations that appear to be reasonable candidates 
for performing numerical evolutions of black hole space- 
times. 



VI. KINEMATIC AL EXTENSION 



Dkij — 



T, = 



+z 7 K ka g a * + z 8 Slg ab K ab , 



(6.2) 



+k igiJ g ab 8 c k + k 5 g k(i 8^g bc + k & g k(l 5 c ])9 ab \ D abc 

+k 7 Tk9i + k 8 9k(iTj), (6.3) 

(6.4) 



kgTi + k w g ab D lab + k n g ab D abi . 



This transformation is (generically) invertable, and is the 
identity when £4 = k\ = kg = 1 and all other za and 
kA vanish. Note that as in KST [l^, it is understood 
that when constructing evolution equations for the trans- 
formed fields u a , all (temporal and spatial) derivatives of 
gij that are introduced by differentiating T a p are to be 
eliminated by substituting the definition of Dkij and the 
evolution equations for g^ . 

We also note that the kinematical transformations 
u a = T a pvP described by Eqs. H6.1|l - (|6.4|) do not change 
the hyperbolicity of the system or the characteristic 
speeds. This is because the characteristic matrix for the 



transformed 



k a 



ZkT a li A k i> v (T- x ) v f i , 



system is £ fc A p = Sfci ^ 
which has the same eigenvalues as £kA p, and the 
the transformed system is S a p = 
' p, which is symmetric and positive 



symmetrizer for 



definite iff is symmetric and positive definite. 

In summary, the unified system of evolution equations 
presented here contains 41 free parameters when written 
in terms of the dynamical fields u a : the 22 parameters 
that entered Eqs. (|3.18() - (|3.24l) as desribed above, plus 19 
transformation parameters {£1, . . . , Zg, k\, . . . , fcn}. In 
Sec. IIVI we reduced the number of free parameters by 6 
to ensure that the system of equations is symmetric hy- 
perbolic. The remaining 35 parameters are freely specifi- 
able and available for other purposes, such as simplifying 
the resulting equations, fixing the characteristic speeds 
to desired values, or optimizing the stability of numeri- 
cal spacetime evolutions. 



Finally we note that the independent dynamical fields 
u a = {g^, Kij, D kij , Q, T. t , N' 1 , M k 1 } can also be 
modified in these evolution equations. It has been 
shown [T2I [l3l l32| that seemingly trivial changes in the 
choice of these dynamical fields can have dramatic ef- 
fects on the stability of numerical spacetime evolutions. 
So following KST [K| we introduce a set of linear trans- 
formations on the dynamical fields. In particular we take 
a new set of fields ii a defined by a transformation of the 
form u a — T a puP, where T a p depends only on various 
parameters and the spatial metric g^ . The most general 
such transformation (which preserves the fundamental 
metric fields {gij, Q, N 1 }) is |51|: 

= + z igi jg ab K ab + z 2 M (i °5 i)a 

+z 3 g t] M a a , (6.1) 
M k l = z A Mk l + z 5 6lM a a + ±z 6 (6t6l-g* a gk b )M a b 
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APPENDIX A: CONSTRAINT EVOLUTION 

The dynamical system derived in Sec. IIIII represents 
the Einstein evolution equations only when certain con- 
straints are satisfied. Here we derive the system of evo- 
lution equations that the constraints themselves satisfy. 
We begin by giving explicit expressions for each of the 
constraints in terms of the dynamical fields and their 
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spatial derivatives: 



c 


= (a ik a jl - a ij a kl )diDi.ii 


+ a i[l a a]h K,4K„ h 




~h\9 9 9 + 2g 9 


.9 + off 1 9 ' 9 






(Al) 


d 


= (g aj g bk St + g c 'g ab 6 k - 


2g ca g^5 k i )D cah K jk 




+2g^H l \d k K ]h 


(A2) 


£i 


= diQ-Tt, 


(A3) 


£k 


= N-'dkN'-Mt 1 , 


(A4) 


£kij 




(A5) 


c 


= 2%^], 


(A6) 


r * 

^nk 


= 2d [n M k] i + 2M [k l (T n] - 


f 2aD n]ab g ab ), (A7) 


Cklij 




(A8) 



We note that these expressions do not contain any spa- 
tial derivatives of the metric fields Q, N a or gif these 
derivatives were replaced by Tj, M k l and Dkij wherever 
they occured. 

The collection of constraints, 

C \C ,Ci, £i, £k , £kij , Cij , C n k ,Cklij}, v^-^) 

is a function of the set of dynamical fields {gij, K-ij, Dkij, 
Q, Ti, N' 1 , Mk 1 }, hence the evolution of the constraints is 
determined by the dynamical evolution equations l|3.18[l - 
(|3.24|) . A straightforward (but very lengthy) calculation 
shows that the evolution of the constraints is determined 
by a system of the form 

d t c a + C ka /3 (u)9 fc c /3 = F a fi{u, du)c p . (A10) 

The coefficients C ka p that appear in this expression de- 
pend on the metric fields Q, N l , and as well as 
the various parameters that enter Eqs. l|3.18[l - i|3.24[) . 
The coefficients F a p depend on all of the dynamical 
fields and their derivatives. The complete expressions 
for the constraint evolution equations are very lengthy, 
so here we list only the principal parts of the equations, 
dtc a ~ —C ka p{u)dkC^ , which can be written: 

d t C ~ N n d n C~^N{2 + 2 X -r 1 )g^d l C J 

+±N{2 + 2^ 10 + 3^ 8 - 4^ 9 ) 5 l ^C, fc fc , (All) 
d t C t ~ -N(l + 2^)d i C + iN(l + 0g jk 9 ab d j C aikb 
-iNg^djCki - iiV(l + 2a)g^g ab d j C kiab 
+N n d n d + ±N(1 - C)9 jk 9 ab d,CkaU, (A12) 
dt£i 0, (A13) 
d t £k l * 0, (A14) 
d t £ kij ~ N n d n £ kij , (A15) 
d t Cij ~ N^d [i C j] +^Nip 2 d k C ij k , (A16) 
d t C nk l ~ N^d 3 C nk l + N^ 7 5\ k d n] C + N^d [n C k]j 
-N^g ab d ln C kh ab ~ N^g ab d [n C k] ab 3 
+N^g ab d [n C Uab \ k] , (A17) 



d t Ci kij ~ 7V a <9 Q C ife „ + Ng a(i d j} Ci k a + A^ 9 <?yd[zCfc]a Q 
+Nx9ijd[fik] - iNri^mdkfij + g^d^i) 
-^N^{g al d {l C k]j a + <7aj9[A]i a ) 
-i7VVio(ft[/^]C ja a + .g^jC^). (A18) 

We note that the principal parts of the constraint evo- 
lution equations can be written in a number of different 
ways using various identities that the constraints satisfy: 



d [i £ j] 

d [k £ 3 ] 1 

d[i£k]ij 
d[kCij] 

d[„C jk ] 1 



d[ n Cki]i 



if ■ 

+M [k l fa + a£ j]ab g ab ) , 

o, 

C[kj % (T n ] + 2aD n ] ab g ab ) 
+M[ n l (Cjk] + 2aC-ik]ab9 ab ) 
+AaM [n l £k ab D 3]ab , 
0. 



(A19) 

(A20) 
(A21) 
(A22) 



(A23) 
(A24) 



We have not yet explored how the character of the con- 
straint evolution equations (e.g. their hyperbolicity) is 
affected by changing their principal parts using these 
identities. 

We note that Eq. i|A10(l implies that the time deriva- 
tive of the constraints vanishes, dtc a = 0, whenever the 
constraints are satisfied, c a = 0, at some initial time. 
This guarantees that the constraints will be satisfied for 
all time in analytic spacetimes. And if these constraint 
evolution equations are symmetric hyperbolic (as we ex- 
pect, but have not yet verified) then this will guarantee 
quite generally that the constraints remain satisfied for 
all time if they are satisfied initially. 



APPENDIX B: THE KST SYSTEM 

The fixed-gauge version of the Einstein evolution equa- 
tions proposed by KST can be shown to be strongly or 
even symmetric hyperbolic for certain choices of the five 
"dynamical" free para meters {a, £, 7, 77, x} character- 
izing that system For the case in which all of the 
adjustable characteristic speeds are set equal to the speed 
of light, we have previously [l3l | determined the regions of 
this parameter space where the equations are symmetric 
hyperbolic, and the regions where they are not. 

Using the technique developed in Sec. IIVI we can now 
construct all of the symmetric hyperbolic fixed-gauge 
KST systems explicitly, even for the general case in which 
the adjustable speeds are left unspecified. For the fixed- 
gauge KST evolution system, Eqs. I|3.1I) - (|3.3I) . the five 
dynamical parameters {a, C, 1, V> x} can be written in 
terms of the symmetrizer parameters {B\, C\, C2, C3, 
C4, C5}. These conditions, which are equivalent to the 
conditions for symmetric hyperbolicity in these systems, 
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ire Specifying the parameters {a, £, 7, 77, %} in this way 

guarantees that the characteristic speeds of this system, 

B 2 = (d+2C 2 )/3, (Bl) 
C = -3C 1 /(C 1 +2C 2 ), (B2) «i = 2a, (B7) 

32d + 10C 2 + 45(C 4 + C 5 + 2gi) «I = - 3C - 4a) - + 6a), (B8) 

7 ~ 135^ ' ( J «§ = 4[(1 + 27)(2 + 2x-tj)-tjC], (B9) 

18d + 45(C 3 + C 4 + 2C 5 ) 

a — — 10(9_B + 2C + 4C ) — ' are real, that the characteristic eigenvectors of the system 

1 are complete, etc. We also note that the expression for 

V = § + B 2 [5(3C 3 + C 4 + 4C 5 ) + 20(C 4 + 3C 5 )a the parameters C and a in Eqs. (jB2j) and pi), enforce 

-3(9C 3 + C 4 + 6C 5 )C] /25(C 3 C 4 - Cf ), (B5) the conditions that C and a be limited to the ranges -3 < 

Q < and < a for the symmetric hyperbolic fixed-gauge 
KST systems no matter what values the characteristic 



X = -§ - J B2[5(C 3 + 2C 4 + 3C5) + 20(2C 4 + C 5 )a 



-3(3C 3 + 2C* 4 + 7C* 5 )C] /25(C 3 C 4 - C 5 2 ). (B6) speeds may have. 
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J. W. York, Jr., in Sources of Gravitational Radiation, 
edited by L. L. Smarr (Cambridge University Press, Cam- 
bridge, England, 1979), pp. 83-126. 

S. Frittelli and O. Reula, Commun. Math. Phys. 166, [22 
221 (1994). 

Y. Choquet-Bruhat and J. W. York, Jr., C. R. Acad. Sci., [23 
Ser. 1 A321, 1089 (1995). 
H. Friedrich, Class. Quantum Grav. 13, 1451 (1996). [24 
S. Frittelli and O. A. Reula, Phys. Rev. Lett. 76, 4667 [25 
(1996). 

C. Bona, J. Masso, E. Seidel, and J. Stela, Phys. Rev. D 
56, 3405 (1997). 

A. Anderson, Y. Choquet-Bruhat, and J. W. York, Jr., 
Topol. Meth. Nonlin. Anal. 10, 353 (1997). [26 
A. Anderson and J. W. York, Jr., Phys. Rev. Lett. 82, 
4384 (1999). [27 
S. Frittelli and O. A. Reula, J. Math. Phys. 40, 5143 
(1999). [28 
M. Alcubierre, B. Briigmann, M. Miller, and W.-M. Suen, 
Phys. Rev. D 60, 064017 (1999). [29 
S. D. Hern, Ph.D. thesis, University of Cambridge (1999), 
gr-qc/0004036. [30 
L. E. Kidder, M. A. Scheel, and S. A. Teukolsky, Phys. 
Rev. D 64, 064017 (2001). [31 
L. Lindblom and M. A. Scheel, Phys. Rev. D 66, 084014 
(2002). 

L. L. Smarr, in Eighth Texas Symposium on Relativis- 
tic Astrophysics, edited by M. D. Papagiannis (The New [32 
York Academy of Sciences, New York, 1977), pp. 569- 
604. [33 
L. I. Petrich, S. L. Shapiro, and S. A. Teukolsky, Phys. 
Rev. D 31, 2459 (1985). [34 
M. W. Choptuik, E. W. Hirschmann, S. L. Liebling, and [35 
F. Pretorius (2003), gr-qc/0301006. 

A. E. Fischer and J. E. Marsden, in General Relativity: 
An Einstein Centenary Survey, edited by S. W. Hawking [36 
and W. Israel (Cambridge University Press, Cambridge, [37 
England, 1979), pp. 138-211. 
L. L. Smarr and J. W. York, Jr., Phys. Rev. D 17, 2529 [38 
(1978). 

M. Alcubierre, Class. Quantum Grav. 20, 607 (2003). [39 
M. Alcubierre, B. Briigmann, D. Pollney, E. Seidel, and 
R. Takahashi, Phys. Rev. D 64, 061501 (2001). [40 
M. Alcubierre, B. Briigmann, P. Diener, M. Koppitz, 



D. Pollney, E. Seidel, and R. Takahashi, Gauge condi- 
tions for long-term numerical black hole evolutions with- 
out excision, gr-qc/0206072 (2002). 

H.-J. Yo, T. W. Baumgarte, and S. L. Shapiro, Phys. 
Rev. D 66, 084026 (2002). 

D. Shoemaker, K. Smith, E. Schnetter, D. Fiske, and 
P. Laguna (2002), in preparation. 

K. AM, Class. Quantum Grav. 19, 5153 (2002). 

C. Bona, J. Stela, J. Masso, and E. Seidel, in Proceedings 

of the Seventh Marcel Grossmann Meeting on General 

Relativity, edited by R. T. Jantzen, G. M. Keiser, and 

R. Ruffini (World Scientific, Singapore, 1996), pp. 410- 

412. 

C. Bona, J. Masso, E. Seidel, and J. Stela, Phys. Rev. 
Lett. 75, 600 (1995). 

O. Sarbach and M. Tiglio, Phys. Rev. D 66, 064023 
(2002). 

M. Alcubierre and B. Briigmann, Phys. Rev. D 63, 
104006 (2001). 

E. Seidel and W.-M. Suen, Phys. Rev. Lett. 69, 1845 
(1992). 

G. B. Cook, M. F. Huq, S. A. Klasky, M. A. Scheel, et al., 
Phys. Rev. Lett. 80, 2512 (1998). 

S. Brandt, R. Correll, R. Gomez, M. Huq, P. Laguna, 
L. Lehner, P. Marronetti, R. A. Matzner, D. Neilsen, 
J. P. E. Schnetter, et al., Phys. Rev. Lett. 85, 5496 
(2000). 

M. A. Scheel, L. E. Kidder, L. Lindblom, H. P. Pfeiffer, 
and S. A. Teukolsky, Phys. Rev. D 66, 124005 (2002). 

D. Shoemaker, K. Smith, U. Sperhake, P. Laguna, 

E. Schnetter, and D. Fiske (2003), gr-qc/0301111. 
G. B. Cook, Phys. Rev. D 65, 084003 (2002). 

J. Balakrishna, G. Daues, E. Seidel, W.-M. Suen, M. To- 
bias, and E. Wang, Class. Quantum Grav. pp. L135-L142 
(1996). 

M. Hoist, to be published. 

M. Shibata and T. Nakamura, Phys. Rev. D 52, 5428 
(1995). 

T. W. Baumgarte and S. L. Shapiro, Phys. Rev. D 59, 
024007 (1999). 

R. Courant and D. Hilbert, Methods of Mathematical 
Physics, Vol. II (John Wiley and Sons, 1962). 
W. H. Press, S. A. Teukolsky, W. T. Wetterling, and B. P. 
Flannery, Numerical Recipes in C (Cambridge University 



11 



[41] 

[42] 
[43] 



[44] 



[45] 
[46] 

[47] 



Press, Cambridge, England, 1992), 2nd ed. 
Y. Choquet-Bruhat and J. W. York, Jr., Banach Center 
Publications 41, 119 (1997). 
T. P. Liu, J. Diff. Equations 33, 92 (1979). 
The factors of the lapse N that appear in this and sub- 
sequent equations were chosen to make these wave equa- 
tions more covariant. Thus for example by using the co- 
efficient nN 2 instead of fi in Eq. 12.31 the characteristic 
speeds relative to a surface normal observer are i^/Jl in- 
stead of ±^fJl/N. 

Note that we can change the principal parts of Eqs. 113.611 
and 13.71 by adding terms proportional to N k (d k Q — Ij.) 
and ^(pkN* - NM k i ). However, doing this for the shift 
equation spoils linear degeneracy, making the system 
genuinely nonlinear. We chose to keep the system lin- 
early degenerate since this condition is known to prevent 
the formation of shocks in similar (but not identical) hy- 
perbolic systems |4^|. 

We note that the 3-index constraint C nk used in this 
paper is not the same as the 3-index constraint Ckij in- 
troduced by KST HI. 

We note that the constraints are to be thought of as 
functions of Q, N l , gij, Kij, Dkij, M k and the first 
derivatives d k Kij, diDkij, dkTi, and diMu ■ The explicit 
expressions for the constraints in terms of these quanti- 
ties are given in Appendix 1X1 

We have considered the most general symmetrizer that 
can be constructed from the spatial metric and constant 
parameters. The symmetrizer given in Eq. 14.21 how- 
ever lacks a number of cross terms between the met- 
ric variables {gij,Q, N z } and the derivative variables 
{Kij, Dkij ,Ti, Mk 1 } ■ These cross terms must vanish as 
a consequence of the symmetry conditions, and so for 
simplicity we omit them from Eq. 14.21 . 



[48] We note that Eqs. KTSl - KTZi are ordered so that each 
of the parameters {B 2 , ipi, . . . , ipw, a, 7, 77, x, C> 
MSi £ ii £ s} appears on the right-hand side of only those 
equations that follow the one in which it appears on the 
left. 

[49] The characteristic speeds associated with these fields are 
actually quite arbitrary. For example we could change 
them to arbitrary multiples of —£, k N /N (including zero) 
simply by adding the terms V k {d k Q-T k ) and V k (d k N i - 
NM k *) to Eqs. |TH§1 and (gjjg) with V k an arbitrary 
multiple of N . Terms of this type do not change the 
symmetric hyperbolicity of the system, but they can af- 
fect its linear degeneracy (and hence the ability of the 
system to form shocks). The vector V k that appears in 
these extra terms could also be specified as an arbitrary 
function of the coordinates. In this case these charac- 
teristic speeds become —£, k (N k — V )/N which can be 
adjusted to arbitrary values with appropriate choice of 
V while leaving the system linearly degenerate. These 
speeds can all be made less than the speed of light in this 
way if that is desired. 

[50] The characteristic speeds associated with the fields Q and 
N l are —£iN z /N and cannot be changed by adjusting the 
various parameters of the theory. These speeds will less 
than the speed of light except when the vector dt becomes 
spacelike. However as noted above the equations can be 
easily modified to make these speeds causal if that is 
desired. 

[51] This is an extension (with a slight change in notation) 
of the transformation introduced by KST [T^. The re- 
lationship between the original KST notation and that 
used here is: z\ = z, k% = k, fca = e, k% — a, ^4 = b, 
k$ — c, and k(, — d. 



